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The ammonia method, recently proposed by Flambaum and Kozlov (2007) to probe variations of the electron- 
to-proton mass ratio, fj, = m c /m p , is applied for the first time to dense prestellar molecular clouds in the Milky 
Way, allowing to test A\i/ \x at different galactocentric distances. High quality radio-astronomical observations 
are used to check the presence of possible relative radial velocity offsets between the inversion transition of NH3 
(J, K) — (1, 1), and the CCS Jn = 2i — lo and N2H + J = 1 — rotational transitions. Carefully selected sample 
of 21 NH3/CCS pairs observed in the Perseus molecular cloud provide the offset AVccs— NH3 — 36 zb 7 s tat i 13.5 S ys 
m s _1 . A similar offset of AV = 40.8 ± 12.9 8tat m s _1 between NH 3 (J, K) = (1, 1) and N 2 H+ J = 1 - has 
been found in an isolated dense core L183 by Pagani et al. (2009). 

Overall these observations provide a safe bound of a maximum offset between ammonia and the other molecules 
at the level of AV < 100 m s _1 . Being interpreted in terms of Afj,/fi, this bound corresponds to Afj,//j, < 1 x 10~ 7 , 
which is an order of magnitude more sensitive than available extragalactic constraints. Taken at face value the 
measured AV shows positive shifts between the line centers of NH3 and these two other molecules and suggest 
a real offset, which would imply a A/i//i~ 4 x 10~ 8 . If A/i/fi follows the gradient of the local gravitational 
potential, then the obtained results are in conflict with laboratory atomic clock experiments in the solar system 
by ~ 5 orders of magnitude, thus requiring a chameleon-type scalar field model. New measurements involving 
other molecules and a wider range of objects along with verification of molecular rest frequencies are currently 
planned to confirm these first indications. 



1. Introduction 

The late acceleration in the universal expan- 
sion reveals a negative-pressure component of the 
bulk of energy density, 'dark energy'. This can be 
a cosmological constant or a dynamically evolv- 
ing scalar field (p, 'quintessence' , which does not 
have a cosmic coincidence problem. However, in 
quintessence the coupling of scalar field to or- 
dinary matter leads unavoidably to long-range 
forces with a variability of the physical constants 
and a violation of the Equivalence Principle. Ef- 
fects that are currently searched in laboratory, 
space-based experiments and astronomical obser- 



vations. 

The suitable constants are the fine-structure 
constant a — e 2 /(Tic) and the electron-to-proton 
mass ratio /1 = m c /m pi or different combinations 
of them with the proton gyromagnetic ratio g p 

a- 

Accurate laboratory atomic clocks experiments 
constraint the temporal variation of a of a/ a = 
(-1.6 ± 2.3) x 10" 17 yr" 1 [19]. Analysis of 
quasar absorption- line systems provided Aa/a = 
(-4.6 ± 1.1) x 10" 6 QI] or |Aa/a| < 6 x 10~ 6 
|22lllj . Here Aa/a = (a' - a) /a, with a de- 
noting the value of the fine-structure constant 
in the laboratory and a' the specific absorp- 



1 



2 



■ 



Velocity dispersion, a v (NH 3 ), km s 



n = 34 


AV = V(CCS) - V(NH 3 ) 


", 1. i 'i, 


V \ 'i 1 


^ 11 1 

f 


n = 21 


1 1 1 1 1 1 1 

1 


: i' , i 'l 


1 \ '■ 1 ' 


j j 

f 



5 6 7 8 

Source velocity, V^, km s~ 



Figure 1. Upper panel: CCS (2i — 1 ) ver- 
sus NH3(1,1) linewidths for cores in the Perseus 
molecular cloud from 20]. The error bars rep- 
resent ler standard deviations. Lower panel: 
selected sample with single-component profiles. 
The dashed and dotted lines are the boundaries 
for pure turbulent and pure thermal line broad- 
ening. 



Figure 2. Upper panel: Velocity offset 

AVccs-NH;, versus the radial velocity for points 
shown in Fig. [I] Lower panel: Same as the upper 
panel but for the points with consistent broaden- 
ing between the lines in Fig 1. The vertical error 
bars include both random and systematic errors 
caused by the uncertainties of the adopted rest 
frequencies. 



tion/cmission line system of a galactic or extra- 
galactic object, the same definition is applied to 
A/i//i. Being linearly extrapolated to rcdshift 
z ~ 2, or t ~ 10 10 yr, the laboratory bound 
leads to |Aa/a| < 4 x 1CP 7 which is below the 
astronomical limits or claims for variability. How- 
ever, laboratory experiments and quasar absorp- 
tion spectra probe different time scales and re- 
gions of the universe, and the connection between 
them is somewhat model dependent |12j . 

The direct laboratory estimate of time varia- 
tion of /j, gives fi/fx = (1.6 ±1.7) x 10" 15 yr" 1 0. 
Astrophysical estimates of changes of /i at high 
redshifs are controversial with claims of nonzero 



and zero changes: A/i//i = (—2.4 ± 0.6) x 10~ 5 
was inferred from the analysis of the H2-bearing 
clouds at z — 2.6 and 3.0 [T7], but this value was 
not confirmed later [7]. More stringent bounds 
are obtained at redshifts z — 0.68 and z — 0.89 
from the analysis of radio-frequency transitions 
in NH 3 and other molecules such as CO, HCO + , 
HC 3 N providing A^i/n = (-0.6 ± 1.9) x 10~ 6 , 
Afi/n = (0.08 ± 0.47) x 10" 6 [US], which favors 
a non- variability of /i at a level of^ 1.5 x 10~ 6 . 

In these studies it is implicitly assumed that 
the rate of time variations dominates over pos- 
sible spatial variations. However, if scalar fields 
trace the gravitational field inhomogcneities spa- 
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tial variation could be present as well [3]. The 
circumstance of spatial variations of constants 
are tested by measurements of the atomic clock 
transition in laboratories on Earth orbiting in 
the changing gravitational potential of the Sun. 
The laboratory experiments constrain the cou- 
plings of a and fj, to the gravitational field at 
a level of fc a)/1 ~ 1CT 6 , 10" 5 0. To reconcile 
these results with quintessence, chameleon- type 
models were suggested which allow scalar fields 
to evolve on cosmological time scales today and 
to have simultaneous strong couplings of order 
unity to matter [5]. In these models the mass 
of the scalar field depends on the local matter 
density, which explains why cosmological scalar 
fields, such as quintessence, are not detectable in 
local tests of the Equivalence Principle. On Earth 
in a dense environment the mass of the field can 
be sufficiently large and ^-mediated interactions 
are short-ranged. On the other hand, the cosmo- 
logical matter density is about 10 30 times smaller 
and the mass of the scalar field can be very low, 
of order of Ho , allowing the field to evolve cosmo- 
logically. 

In this paper we apply the ammonia method for 
the first time to Milky Way sources. A detailed 
and exaustive account has been provided in [10] . 

2. Interstellar molecules probe A////Z spa- 
tial variations 

Among molecules, ammonia is of particular in- 
terest due to its high sensitivity to changes in the 
electron-to-proton mass ratio, /i [4]. For NH3, the 
sensitivity coefficient of the inversion transition 
v = 23.69 GHz was calculated in [?]: 

^=4.46^. (1) 

V [1 

By comparing the observed inversion frequency of 
NH3 (1,1) with a suitable rotational frequency of 
another molecule arising co-spatially, a limit on 
the spatial variation of [i can be determined: 

Am = Q 289 F rot -ti nv ^ 28g AF 

fl c c 

where V mt and Vi nv are the apparent radial veloc- 
ities of the rotational and inversion transitions, 
respectively. 



Thus the comparison of the relative radial ve- 
locities of ammonia inversion lines and rotational 
transitions of another N-bearing (N 2 H + ) and C- 
bearing (CCS, HC3N) molecules can be used 
to set a limit on spatial variations of A/if/j,. 
These molecules are often detected in the molec- 
ular clouds allowing to probe this constant in 
the Milky Way. In molecular clouds consider- 
able fraction of material is not participating in 
the process of star formation residing in cold 
and dense clouds which are well suited for ac- 
curate measurements. These can be low-mass 
protostellar cores containing infrared sources, or 
even starless cores without embedded luminous 
sources. The physical conditions of the latter 
are of dense molecular gas preceding gravitational 
collapse and when the gas density in the core cen- 
ter is n H2 < 10 5 cm~ 3 they can be dynamically 
stable against gravitational contraction [5] . Typ- 
ical cores have mean gas density n H ~ (1 — 2) x 10 4 
enr 3 , velocity dispersion a v — 0.17 km s _1 , ki- 
netic temperature Tkin = 11 K, radius R — 0.09 
pc, mass M ~ 1 M Q and NH 3 , N 2 H+, CCS, and 
HC3N are usually observed in emission. 

3. Cores in the Perseus molecular cloud 

An ammonia spectral atlas of 193 dense pro- 
tostellar and prestellar cores of low masses in 
the Perseus molecular cloud is provided by |20j . 
The spectral observations of the cores in NH3 
(1,1), (2,2), CCS (2 X - 1 ) and CC 34 S (2 X - 1 ) 
lines were carried out simultaneously using the 
100-m Green Bank Telescope (GBT). Each tar- 
get was observed in a single-pointing, frequency- 
switched mode. The GBT beam size at 23 GHz 
is FWHM = 31" or 0.04 pc at the distance of 
260 pc of Perseus cloud. The spectral resolution 
was 24 m s _1 and both molecules were observed 
with similar angular and spectral resolutions. 

The atlas of 193 sources is very conveniently 
presented on the website (see [5U] for details) 
where the original spectra are shown along with 
the best fitting models, model parameters and 
their uncertainties. The number of cores where 
both CCS (2! - 1 ) and NH 3 (1,1) lines were 
detected is 98. The central velocities of these 
lines, being averaged with weights inversely pro- 
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portional to the variances of the measurements, 
reveal that AVccs-nh 3 = 16 m s _1 . It was sug- 
gested [5D] that this offset is due to uncertainties 
in the assumed rest frequency of the CCS line 
which is of 13.5 m s _1 (cfr Table 1). However, 
the unweighted averaging, which is more appro- 
priate in the presence of unaccounted errors, gives 
AV n=gs = 44 ± 13 m s" 1 

Many cores show evidence for the multiple ve- 
locity components along the line-of-sight. The 
selection of simple profiles which have been satis- 
factorily fitted by a single component model pro- 
vides a sample of n — 34 pairs The weighted mean 
for this sample gives AV™ =3i = 40 ± 10 m s _1 . 
and AV n= 34 = 39 ± 10 m s _1 for the unweighted 
mean [TO]. 

A further refinement of the sample could be 
done taking only cores which show line broaden- 
ing values between pure turbulent, cr„(CCS) = 
<7„(NH3), and pure thermal, cr„(CCS) = 
0.55cr^(NH 3 ), regimes, which should be the case 
if both emission lines arise co-spatially. This re- 
quirement further reduce the sample to n — 21 
clouds with AV™ =21 = 50 ± 0.013 m s" 1 , and 
AV n= 2i — 48 ± 13 m s _1 for the weighted and 
unweighted mean respectively. 

The fact that scatter remains rather constant 
in the various samples reflects the complexity of 
gas kinematics and effects of chemical segregation 
of one molecule with respect to the other in dense 
molecular cores which are the dominant sources 
of error. 

The unweighted mean from the n — 21 sam- 
ple could be taken as the best estimate of the 
velocity offset for the Perseus dark cores, but 
in fact we could have taken any of the former 
samples without changing significantly the result. 
When we account for uncertainties in the rest fre- 
quencies , we have AVccs-nh 3 =48 ± 13 ± 13.5 
m s _1 . Other observations in the Pipe neb- 
ula JT5] and in Infrared dark clouds |H] are 
examined in [10]. 8 pairs in the Pipe Neb- 
ula provide AVccs-nh 3 =53 ± ll sta t ± 13.5 sys 
m s- 1 . 36 NH 3 /N 2 H+ and 27 NH3/HC3N pairs 
observed in Infrared Dark Clouds, which are 
much more massive than the dense cores, provide 
AVn 2 h+-nh 3 = 148 ± 32 sta t ± 13.6 sys m s" 1 and 
AVHC 3 N-NH 3 =115±37 8ta t±31 sys m s" 1 , respec- 



tively, thus supporting the results found in the 
Perseus region. 

4. The L183 dense core 

Recently Pagani et al. [TB] have compared the 
N2H + and N 2 D + transitions with the transitions 
of NH 3 in the dark cloud LI 83. The whole elon- 
gated dense core of L183 has been fully mapped in 
N 2 H+ and N 2 D+ (J: 1-0) lines with the IRAM 30- 
m telescope with velocity resolution in the range 
30-50 m s _1 . Spatial resolution ranges from 33 
arcsec at 77 GHz to 9 arcsec at 279 GHz. The spa- 
tial sampling is 12 arcsec for the main prestellar 
core and 15 arcsec for the southern and northern 
extension of the prestellar core. Observations of 
NH3 (1,1) and (2,2) inversion lines have been per- 
formed at the Green Bank 100-m telescope with 
velocity sampling of 20 ms _1 , angular resolution 
of ^35 arcsec and spatial sampling of 24 arcsec 
all over the source. 

Extremely narrow molecular lines observed in 
this cold dark cloud provide a sensitive spectro- 
scopic tool to track small systematic velocity gra- 
dients and to align values of poorly known fre- 
quencies with more reliable ones. For the latter, 
NH3 is the molecule of choice since it is stable and 
its frequencies are accurately measured in the lab- 
oratory. NH3 (1,1) average frequency from the 
whole HFS is v = 23 694495487 (±48) Hz [5], 
or v = 23 694495481 ± 22 Hz from the revision 
of Hougen [T3]. Thus NH3 provides a reference 
transition known with a precision of few 10~ 10 or 
~ 0.6 m s _1 . 

Ammonia and diazenylium have the same 
chemical origin, starting from N 2 and are well- 
known to be coexistent as discussed by e.g. 
23 2j]. In L183 the velocity along the dense fila- 
ment is constantly changing suggesting a flow to- 
wards the prestellar cores and a rotation of the fil- 
ament around its vertical axis. NH 3 (1,1), N 2 H + 
and N 2 D + (J: 1-0) all trace exactly the same gra- 
dients and it seems therefore compulsory that 
their velocities be identical as there is no obvious 
possibility that the velocity gradients be exactly 
parallel but offset from each other, especially in 
the case of a cylinder rotation. Finally, the fact 
that N 2 D + velocity centroids are almost identi- 
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Table 1 

Molecular transitions and uncertainties, e v , The 
numbers in parentheses correspond to la errors. 



Transition Vteat > s v , 

GHz m s~ 



CCS J N = 2i - 1 22.344033(1) 13.4 

NH 3 (J, K) — (1, 1) 23.694495481(22) 0.6 

HC 3 N J = 5 — 4 45.4903102(3) 2.8 

N 2 H+ J = 1 -0 93.173777(4) 13.5 



cal with those of N 2 H + indicates that the differ- 
ent opacities of the lines are not introducing any 
measurable bias here. 

A similar behaviour of all the three molecules 
shows that these share the same volume of the 
cloud and undergo the same macroscopic velocity 
shifts. 

On the basis that all the three species are spa- 
tially coexistent and trace the same velocities Pa- 
gani et al [16] used the observed offset to adjust 
the frequencies of N 2 H + and N 2 D + to that of 
NH 3 . For the reference position alone, the differ- 
ence is of AVn 2 h+-nh 3 =40.8 ± 0.56 m s _1 and 
AVn 2 h+-nh 3 =40.8 ± 12.9 m s _1 on 65 common 
positions . For the three transitions of N 2 D + , the 
direct comparison of the reference position with 
NH 3 (1,1) spectrum yields (J:l-0) : 37.7 ± 0.85 
m s- 1 (J:2-l) : 47.7 ± 0.92 m s" 1 (J:3-2) : 63.6 ± 
4.7 m s _1 and similar offsets are found for N 2 D + 
(cfr [16] for details). 

The rest frequencies for the molecule N 2 H + in 
|16j analysis were taken from [T] and are known 
with precision of ~ 7 kHz. However, an error 
of 4 KHz is given in the Cologne Database for 
Molecular Spectroscopy (CDMS). The difference 
between the data from [1] and the CDMS frequen- 
cies is — 12 m s" 1 and the Pagani et al velocity off- 
set between N 2 H + and NH 3 should be corrected 
by this difference when adopting the CDMS rest 
frequency. 

5. Discussion 

5.1. A safe bound for Afi/fi 

At the high spectral resolution of these ob- 
servations the errors in the molecular line posi- 



tion measurements arc mainly restricted by the 
uncertainties in laboratory frequencies, e v ~ 1 
kHz, which correspond to the Vlsr uncertain- 
ties of e v ~ 10 m s _1 Taking into account that 
A/x///~ 0.3AV/c, these errors implies a sensitiv- 
ity in A/Lt//it at a level of ~ 10~ 8 , i.e. about 100 
times more sensitive than the Aji/ \i estimate de- 
duced at z = 0.68 [4]. This level of accuracy can 
be achieved when assuming that molecules are 
co-spatially distributed within the cloud, and are 
observed simultaneously with the same receiver, 
beam size, system temperature, and velocity reso- 
lution. Violation of any of these conditions leads 
to shifts of the line centers. Overall the obser- 
vations in the Perseus dense cores as well as in 
the L183 core allow a maximum offset between 
ammonia and the other molecules of AV < 100 
m s _1 in the most conservative way, which corre- 
sponds to A/i//i < 1 x 10~ 7 , i.e. an order of mag- 
nitude more sensitive than previous astronomical 
constraints on /i. 

This firm bound is in conflict with the value 
Afi/fi = (—24 ± 6) x 10~ 6 obtained from molec- 
ular hydrogen H 2 absorption lines at z ~ 3 |17j . 
However, the variability of fi was not confirmed 
at the level of \Afi/fi\ < 4.9 x 10~ 5 [25] and, re- 
cently, a more stringent limit on A/i//i was found 
at z ~ 3, Afx/fj, = (2.6 ± 3.0) x 10" 6 [7]. 

These extragalactic measurements of Afi/fj, are 
performed in systems in which gas densities are 
similar to those in the Milky Way clouds. Thus 
also in the high-z absorbers the dependence of \i 
on p is extremely weak, and the value of A/1//1 
in quasar absorbers is expected to be at the same 
level as in the interstellar clouds, i.e. ~ 10~ 7 . 
If a Afi/fi is found above this value, then the 
variation should be temporal rather than spatial. 

5.2. Hints of variation of A/i//i ? 

Taken at face value the measured AV^ in 
Perseus and L183 being similar while involving 
different molecules show the possibility of a pos- 
itive shifts between the line centers of NH 3 and 
other molecules of the order of 40 m s _1 . This 
would imply a A/x///~ 4 x 10~ 8 . 

In theories which consider a direct coupling 
between the local matter density and the scalar 
fields driving varying constants the spatial change 
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of the constants is proportional to the local 
gravitational potential [15] (and references cited 
therein). The effect of scalar fields on p can be 
expressed as 

^ = k,A$ , (3) 
A* 

where k^ is a dimensionless coupling constant of 
a very light scalar field to the local gravitational 
potential <I>, and A<i> is the difference of the grav- 
itational potentials between two measurement 
points. Since in our case A<i> ~ 10~ 7 , a value 
of Afi/ fi= 4 x 1CP 8 would require k^ ~ 1, which 
is inconsistent with the fc M < 10 _5 estimated from 
atomic clock experiments [5] . 

However, chameleon models are specifically 
conceived to circumvent this type of conflict. 
The baryon masses and coupling constants are 
strongly dependent on the local matter density, 
p, [T5]- The dynamics of the scalar field in these 
models depends on p: in low-density environ- 
ments it is determined by the scalar field po- 
tential V(ip), whereas in high-density regions like 
the Earth's surface it is set by the matter-^ cou- 
pling. The range of the scalar-mediated force 
for the terrestrial matter densities is then very 
short, less than 1 mm [IB] . Because of the strong 
(^-coupling to matter the measurements of the 
frequency drifts in the atomic clock experiments 
ought to be insensitive to the changes in the grav- 
itational potential at Earth caused by the eccen- 
tricity of Earth's orbit. 

We note that the difference of matter densi- 
ties between the terrestrial environment (p^ ~ 
3 x 10 24 GeV cm~ 3 ) and, e.g., dense molecu- 
lar cores (pdoud ~ 3 x 10 5 GeV cm" 3 ) is of 
20 orders of magnitude. Thus, Ap/p variation 
Ap/p = 4 x fCP 8 require models which treat the 
(^-mediated force as a short-range force depend- 
ing on the matter density. In accord with model 
predictions, the estimated value of the scalar field 
is much less than one. In this way the measured 
spatial variation of Ap/p does not contradict the 
laboratory studies on atomic clocks due to the 
extremely different density environment in the 
terrestrial measurements and in the interstellar 
medium. 



6. Conclusions 

1. Several molecules and in different envi- 
ronments in the Milky Way are used to 
probe Ap / p variations by means of the 
NH3 method. This is done by using CCS 
measured in almost one hundred of dense 
Perseus clouds and N2H" 1 " and N2D + in a 
detailed high spatial resolution survey of 
the dense core LI 83. 

2. These observations provide evidence that 
NH3 does not change by more than 100 
m s _1 allowing us to place a robust bound 
of Ap/p< 10- 7 which is more than an or- 
der of magnitude more stringent than the 
extragalactic bounds of the same quantity. 

3. Both sets of comparison point for a com- 
mon offset of about 40 m s _1 . Although 
it is possible that this offset results from 
rest frequency uncertainties the fact that 
several transitions are involved suggest that 
the effect might be real. An intrinsic shift 
of 40 m s _1 would correspond to Ap/p~ 
4 x 10~ 8 . 

4. To cope with solar system experiments the 
effect would require chameleon-type scalar 
field models which predict a strong depen- 
dence of masses and coupling constants on 
the ambient matter density. 

5. New observations to verify this possibility 
are solicited, and some are underway at ra- 
dio telescopes of Medicina, Effelsberg and 
Nobeyama. 
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